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We compute the differential distributions of two three-jet event-shape observables 



^ I at the next-to-leading order accuracy at fixed values of the DIS kinematic variables. 

' The observable i^out measures the out-of-event-plane momentum. The other observ- 

^.f^ . able ?/3 is the maximum value of the ycut resolution variable for which an event is 

classified as three-jet event. We also show the dependence of the fixed-order predic- 
. ^ tions on the renormalization and factorization scales. The radiative corrections are 

. in general large and depend on the value of the DIS kinematic range. 
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^ ■ 1 Introduction 

The analysis of event-shape observables in e"^e~-annihilation and in deeply- 
inelastic lepton-proton scattering (DIS) proved to be a powerful method to 
study Quantum Chromodynamics (QCD) [1]. The standard QCD analysis 
of event-shapes consists of matching the next-to-leading order (NLO) and 
resummed next-to-leading logarithmic (NLL) predictions, possibly improved 
with analytic predictions for the power corrections (PC). Such matched pre- 
dictions also describe the distributions of multi-jet rates in e~^e~ annihilation 
with high accuracy [2]. Thus, it is interesting to investigate whether similar 
level of accuracy can also be achieved in predicting distributions of multi-jet 
observables in DIS. 

Two-jet ^ event-shape observables in DIS have been thoroughly analyzed both 



^ In counting the number of jets in this paper, we neglect the ubiquitous beam-jet. 
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theoretically [3] and experimentally [4]. Three-jet rates have been also ana- 
lyzed [5] based on computations valid at the NLO accuracy [6]. Rcsummed pre- 
diction to the three-jet rates have not yet been computed, therefore matched 
predictions are not available for jet rates. However, much progress has been 
achieved in computing resummed predictions at the NLL accuracy for three- 
jet event shapes [7,8,9] that arc sensitive to large angle soft emission and thus 
exhibit rich geometry- dependent structure. These developments in theory en- 
couraged the experimenters to consider three-jet event shapes [10]. However, 
predictions to three-jet event shapes at the NLO accuracy have not been com- 
puted yet. In this letter we aim to fill this gap. 

It is well-known that the validity of fixed-order predictions is rather con- 
strained. In the perturbative expansion of the distribution of an observable 
O, logarithmic terms log" O appear which require all-order resummation if 
the value of O is small. In analyzing the data, most of the statistics he in the 
range of small O, thus the computation of resummed predictions is indispens- 
able for the experimental analysis. The observables we choose to compute are 
those for which resummed predictions are known. 

The first three-jet event-shape observable in DIS that has been computed at 
the NLL accuracy is the /Cout variable, that measures the out-of-event-plane 
QCD radiation. The observable i^out was defined in Ref. [7] as the sum of the 
momentum components perpendicular to the event plane, 

^out = Ebri- (1) 

h 



The summation extends over all particles (hadrons in the experiment, partons 
in the theoretical computation). The event plane is spanned by the proton 
three-momentum p and the unit vector n that defines the thrust major axis 
in the plane perpendicular to the beam, 

Tu = max j{^\'Ph-n\, n-p^O. (2) 



Recently the caesar program has been published [8,9] that can be used for 
computing cross sections of two- and three-jet event shapes ^ in an automatic 
way. In particular, the distribution of the ys observable, that is defined to 
be the largest value of the jet resolution variable t/cut such that the event is 
clustered into three jets, is also known to NLL accuracy [9]. For defining the 
jets, the computation uses the A;_L-clustering algorithm of Ref. [11]. 

In this paper we compute the distributions of the two three-jet event shape 
^ It can also be used for computing dijet event shapes in hadronic collisions. 
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observablcs in DIS for which resummcd predictions already exist, but the 
fixed-order radiative corrections have not been computed before. We use the 
NLOJET-)--)- program [12]. 



2 Computation of fixed-order predictions 



The next-to-leading order cross section for electron-proton scattering into 
three jets is the convolution of the parton density function of the incoming 
proton and the hard scattering cross section, 



(3) 



where p'^ and are the four-momenta of the incoming proton and the ex- 
changed virtual photon, respectively, fa/piv^t^'F') is the density of the parton 
of type a in the incoming proton at momentum fraction rj and factorization 
scale /If- The corresponding parton level cross sections are 

q)^l da^ip, q) = J dr(3)(|Mf ^ J^^) , (4) 

3 3 



and the next-to-leading order correction is sum of three terms 



= J da^ip, q)+ j daj(p, q) + j d(7„C(p, q) , (5) 

4 3 3 

where dcr^ and dcr^ are the real and virtual contributions to the partonic cross 
section. The contribution da^ represents the coUinear-subtraction counter 
term. The pole structure of this term is unique, while its finite part depends 

on the factorization scheme. We use the MS scheme as defined precisely in 
Ref. [13]. The parton density functions are also scheme dependent, so that the 
scheme-dependence cancels in the hadronic cross section of Eq. (3). 

The three integrals on the right hand side of Eq. (5) are separately divergent 
but their sum is finite provided the jet function J^"*^ defines a coUinear and 
infrared safe quantity, which formally means that 

J(4) _^ J{3) ^ (6) 
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whenever the four-parton and threc-parton configurations arc kincmatically 
degenerate (regions of one unresolved parton). In addition, to define a three- 
jet observable, both J^^^ and J^^^ have to vanish if less than three partons are 
resolved. The presence of the singularities means that the separate pieces have 
to be regularized and the divergences have to be cancelled. We use dimensional 
regularization m d = A — 2e dimensions in which case the divergences are 
replaced by double and single poles of the form and 1/e. We assume that 
ultraviolet renormalization of all Green functions to one-loop order has been 
carried out, so the divergences are of infrared origin. In order to get the finite 
sum one has to rearrange the various contributions by subtracting and adding 
the same, in = 4 dimensions singular terms to the three contributions in 
Eq. (5) so that each becomes separately finite in o? = 4 dimensions. 

The essence of this rearrangement is to define a single subtraction term du^ 
that regularizes the divergences in the real term which comes form the unre- 
solved soft and collinear regions. Thus, the three singular integrals in Eq. (5) 
are substituted by three finite ones: 



where the four-parton integral is given by 

= / [d<7„^(p, g).=o - datip, g).=o] . (8) 



We have two three-parton contributions to the NLO correction. The second 
term on the right hand side of Eq. (7) is the sum of the one-loop contribution 
and a Born term convoluted with a universal singular factor /, 



The factor / contains all the e poles which come from the dcr^ and do"^ 
terms that are necessary to cancel the (equal and with opposite sign) poles 
in dcr^. The ® operation means correlations in colour space. The last term in 
Eq. (7) is a finite remainder, in the form of a convolution, that is left after 
factorization of initial-state collinear singularities into the non-perturbative 
parton distribution functions. 



1 

/da:ar°W(x,xp,g) 
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1 

5: / dx I [da^,{xp, q) ^ [P{x) + K{x)r\^^ , (10) 

«' 3 

where the x-depcndcnt functions P and K are similar (but finite for e ^ 0) 
to the factor I. These functions are universal, that is, they are independent 
of the scattering process and of the jet observables. 

There are many ways to define the da^ subtraction term, but all must lead 
to the same finite next-to-leading order correction. In computing the NLO 
corrections to multijet cross sections the dipole subtraction scheme of Catani 
and Seymour [13] is a convenient formalism. It is used both in the disent and 
the NLOJET-I— I- programs. The subtraction scheme applied in the nlojet-|— |- 
program is modified slightly as compared to the original one in [13] in order 
to have a better control on the numerical computation. The main idea is to 
cut the phase space of the dipole subtraction terms as introduced in Ref. [14]. 
We thus define the da^ local counter term as 



da^ = ^dr(^)(p„,g,Pi,...,P4) ^ 



{4} -^{4} 
H H Dij,k{Pa, q,Pl, ■ ■ ■,Pi)J^^\Pa, ■ ■ ■,Pij,Pk, ■ ■ ■)Q{yij,k < a) 

L pairs k^i,j 
i,j 

+ Yl ^tjiPa^ (i^Pi^ ■ ■ -^PdJ^^^Pa, ■ ■ ■ ,Pij, ■ ■ O©!! - ^ij,a < «) 

pairs 

+ T.Dt(Pa, q,Pi,.. .,pdJ^''\Pa, ...,Pk,.. -Mui <a)\, (11) 

i^k ) 

where dF^^^ is the four-parton phase space including all the factors that are 
QCD independent, Z^{4} denotes the sum over all configurations with 4 partons 
and S'{4} is the Bose symmetry factor of the identical partons in the final state. 
The Dij^ki and functions are the dipole factors given in Ref. [13]. The 
function J^^^ is the jet function which acts over the three-parton dipole phase 
space. The variables Uij^k, ^ij,a and Ui are the dipole variables used for defining 
the exact factorization of the phase space [13]. The parameter a e (0,1] 
controls the volume of the dipole phase space. The case of a = 1 means the 
full dipole subtraction. We can speed up the computer program by choosing 
a ~ 0.1, which keeps the subtraction in the vicinity of the singular regions, but 
avoids the CPU- intensive computations of the dipole terms where those are not 
necessary. Furthermore, checking that the predictions are independent of the 
parameter a, that sets the volume of the cut dipole phase space, gives a strong 
control that indeed the same quantity is subtracted from the real correction 
as added to the virtual one. Choosing a < 1, the insertion operators /(q;,^), 
P(q;), K{a) depend on a. The explicit expressions can be found in Ref. [15]. 
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The integrand of the NLO contribution (7^'"0{4} ^-g^ contains integrable 

square root singularities. Integrating these singularities by simple Monte Carlo 
integration technique (choosing random values of the integration variables 
uniformly) is not efficient because the variance of the estimate of the integral 
is formally infinite, therefore, the estimate of the statistical error of the integral 
is unreliable. To improve the convergence of the Monte Carlo integral, in the 
NLOJET++ program the phase space is generated by multi-channel weighted 
phase space generator [16]. 

Once the phase space integrations are carried out, we write the NLO jet cross 
section in the following form: 

1 

'^^■^^ = / d?7 fa/p(V, /^f) C^SlO (Pa, Q, Ois{^l%) , ^4t/ QhS, ^^F/QHs) , (12) 

where (T^'^lq represents the sum of the LO and NLO contributions to the 
partonic cross section, given in Eqs. (4) and (7) respectively, with jet function 
J. In addition to the parton momenta and possible parameters of the jet 
function, it also depends explicitly on the renormalized strong coupling q;s(a*r), 
the rcnormalization and factorization scales /x^ = x^Qh.s. ^-^d = xfQu.s., 
where Qu.s. is the hard scale that characterizes the parton scattering. The scale 
Qh.s. is usually set event by event to a measurable energy scale of the event. 
Furthermore, the cross section also depends on the electromagnetic coupling, 
for which we used MS running q;bm(Q^) at the scale of the virtual photon 
momentum squared, — —(f. 

The publicly available version of the nlojet-| — \- program [12] is based on the 
tree-level and one- loop matrix elements given in Refs. [14,17], crossed into the 
photon-parton channel. It uses a C/C-l— I- implementation of the LHAPDF 
library [18] with CTEQ6 [19] parton distribution functions and with the cor- 
responding Og expression for the renormalized coupling which is included in 
this library. The CTEQ6 set was fitted using the two-loop running coupling 
with a^{Mzo) = 0.118. 



3 Results 

We computed the distributions for fixed values of the DIS kinematic variables 
and xb, as done in the resummation computations [9]. We used three sets 
of values: = (35 GeV)^ with xb = 0.02 and 0.2 and = (65 GeV)^ with 
Xb = 0.2. The = (65 GeV)^ with xb = 0.02 kinematical point is not 
accessible at HERA. In order to select events with at least two jets with high 
transverse momenta {p± ~ Q), we required y2 > 0.1, where y2 is the maximum 
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value of the jet resolution variable Ucnt such that the event can be classified as 
a two-jet one [7]. We employed a cut on the rapidity of the final-state momenta 
in the Breit frame such that 



Figs. 1 and 2 show the differential distributions for -f^out and observables. 
The shaded bands correspond to the range of scales 1/2 < < 2, where 
~ I^r/Qh.s. ~ A*f/Qh.s.' with hard scattering scale chosen to be Qh.s. ~ Q^- 
We observe several features of these plots. 

Let us consider first the distributions in A'out- We find that the radiative cor- 
rections are in general large, thus the scale-dependence reduces only relatively 
to the cross sections. The corrections also depend strongly on the values of the 
DIS kinematic variables: they increase with decreasing and with decreasing 
xb (as seen from the plots in the right panel). They also increase with increas- 
ing value of A'out because the phase space for events with large out-of-plane 
radiation with three partons in the final state (at LO) is much smaller than 
that with four partons in the final state (real corrections). The boundary of 
the phase space in Kout depends on the value of xb, increases with decreasing 
Xb, and is about 20% larger for the NLO computation than at LO. The cross 
sections decrease rapidly with increasing /Tout- The rate of this decrease also 
depends on xb, being much quicker for larger values of Xb due to the smaller 
phase space. Nevertheless, the small cross section for medium or large values 
of Kout leaves the small Xout-region for experimental analysis. 

In the small i^^out-region, the logarithmic contributions of the type \nKont/Q 
are dominant as can be seen on the plots in the right panel. At LO, the 
cross section behaves as —as{Q) InKout/Q, while at NLO the asymptotically 
leading term is oi'^{Q) In^ K^ut/Q for small values of K^ut, therefore, the fixed- 
order predictions diverge with A'out with alternating signs, which makes 
the resummation of these large logarithmic contributions mandatory. Reliable 
theoretical predictions can be obtained by matching the cross sections valid 
at the NLO and NLL accuracy as described in Ref. [20,3]. 

Similar qualitative conclusions can be drawn from the distributions with 
some important differences. Although the corrections are also large, they are 
much smaller than in the case of the ATout- Thus the reduction in the scale- 
dependence can clearly be seen. An important reason for the smaller correc- 
tions is that the phase space in is the same at LO as at NLO. The size of 
the corrections depends more on the value of the momentum transfer, less on 
Xb- However, the phase space in y^ is more dependent on xb (increases with 
decreasing xb), than on Q"^ (decreases slightly with increasing Q"^). The cross 
section is sizable again for small values of the event shape and the need for 




(13) 
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Xb = 0.2, Q = 35 GeV 
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Fig. 1. The differential distribution of the i^out observable at three different fixed val- 
ues of the DIS kinematic variables. The left panel shows the distributions as a func- 
tion of Ko^t/Q, the right panel shows the distributions as a function of Ini^'out/Q in 
order to exhibit the logarithmic dominance for small values of the observable. The 
LO predictions are shown with dashed lines, the NLO predictions are shown with 
solid lines. The errorbars indicate the uncertainty of the numerical integration that 
is negligible for the computation at LO accuracy. 
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resummation in this region is clearly seen on the right panel. 



4 Conclusions 

In this paper we presented a computation of NLO corrections to the differential 
distributions of the three-jet event-shape observables i^out and in DIS. We 
found large radiative corrections, especially in the case of -R'out, indicating that 
the inclusion of even higher order corrections as well as the non-perturbative 
power corrections is necessary in order to make a reliable prediction, useful 
for experimental analysis. The cross sections decrease rapidly with increasing 
values of the event-shape variables, leaving the region of small values of the 
observables with sufficient statistics for an experimental analysis of data col- 
lected at HERA. In these regions the all-order resummed predictions in the 
NLL approximation are known, therefore, the matching of the NLO and NLL 
distributions promises us reliable predictions. Such results are expected to be 
available soon [21]. 

We are grateful to G. Zanderighi and A. Banfi for their helpful correspondence 
on DIS event-shape observables. This work was supported by the Hungarian 
Scientific Research Fund grant OTKA T-038240 and by the Swiss National 
Science Foundation (SNF) under contract number 200020-109162. 
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